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Abstract. Travelling modes of single-photon-added coherent states (SPACS) are 
characterized via optical homodyne tomography. Given a set of experimentally mea- 
sured quadrature distributions, we estimate parameters of the state and also extract 
information about the detector efficiency. The method used is a minimal distance 
estimation between theoretical and experimental quantities, which additionally allows 
to evaluate the precision of estimated parameters. Given experimental data, we also 
estimate the lower and upper bounds on fidelity. The results are believed to encourage 
preciser engineering and detection of SPACS. 

PACS numbers: 03.65.Ta, 03.65.Wj, 42.50.Xa, 42.50.Dv 
1. Introduction 



Optical homodyne tomography is a powerful technique to infer continuous-variable 
quantum states of a specific mode of electromagnetic radiation. Its history saw a 
dramatic boom in the last two decades, when both theoretical and experimental 
methods evolved significantly from the first proof-of-principle studies [U [2j [3] to the 
state-of-the-art detection of arbitrarily shaped ultrashort quantum light states [I] and 
the experimental analysis of decoherence in continuous- variable bipartite systems [5]. 
Different stages of the research in this area can be seen in the monographs and 
reviews [6j [TJ [9], [101 [11], where different approaches to the state reconstruction are 
outlined and corresponding experimental realizations are discussed. 

The goal of this paper is to consider both theoretically and experimentally the 
detection of single-photon-added coherent states (SPACS). These states are defined 
by formula a*\oi) / ^1 + |a| 2 , where \a) is a conventional coherent state (a G C) 
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and is a photon creation operator. Photon-added states of light and their 
non-classical properties were considered originally in the papers [T2l [13] and then 
realized in practice [HI EE]. The techniques of photon addition and subtraction 
allowed to check experimentally the commutation relation between the corresponding 
operators [T6J [T7J HE]. The process tomography of photon creation and annihilation 
operators was recently reported [19]. Nonclassical behaviour of photon-added states 
was demonstrated in the papers [201 12H 122] and noiseless amplification was discussed in 
Ref. [23]. 

The practical homodyne detection of some signal results in experimental quadrature 
distributions w ex (X, 9) to be compared with theoretically predicted ones Wth(X, 9). The 
adequate theoretical model should take losses into account, which are usually modelled 
by fictitious beamsplitters with transmittivity rj placed in front of ideal detectors. In 
the paper [15], the explicit form of such theoretical quadrature distributions u> t h(A, 9) 
for SPACS is found for any 77. We can associate density operators p ex and p t h with 
distributions w ex (X, 9) and w t h(X,9), respectively. Note that these states are mixed in 
general and depend on parameters a and rj. 

We can naturally define the fidelity of detection as fidelity between p ex and p t h, i-e. 
F = Tr| v /p^y / /^h| = Tr *J y/p~thp cx y/pih- It is tempting to express F directly through the 
measured distributions w ex (X, 9) avoiding reconstruction of the state p ex and dealing 
with complicated formulas. The easiest way is to find the Bhattacharyya coefficient [21] 
for distributions w ex (X, 9) and w t h(X, 9), which turns out to be the upper bound for 
F [25]. Alternatively, one can use upper and lower bounds for F 2 developed in the 
paper [26] and also known as super- and sub-fidelity, respectively. In this paper, we 
present operational ways to calculate these quantities. 

In principle, maximizing the sub-fidelity with respect to a and 77 would enable us 
to estimate both these parameters. As it will be shown by an example in Section HI 
such a method can be applied to extremely precise data only. If this is not the case, 
parameters a and i] can be estimated by minimizing another distance between the states 
p ex and pth (not the Bures distance related to the fidelity). Fortunately, the Hilbert- 
Schmidt distance is easy to compute via tomograms and its minimization is performed in 
Section [3j As a result, an operational estimation of state and measurement parameters 
is achieved. Finally, errors of the estimated parameters are evaluated by using symmetry 
condition w(X, 9 + tt) = w(—X, 9) met by fair optical tomograms. This approach was 
suggested and demonstrated in the papers [27] and [2E], respectively. The improved 
precision of homodyne detection is of vital importance to check different uncertainty 
relations (see [28] and references therein) as well as to probe commutation relations 
between position and momentum of massive particles, which may be modified by gravity 
and feasibly detected with the help of quantum optics [29] . 

The paper is organized as follows. 

In Section |2j we remind the explicit formula of homodyne quadrature distributions 
of SPACS modified by the losses. In Section El we present theoretical basics and 
demonstrate particular results of the minimal distance estimation of the state and 
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apparatus parameters. In Section [U the fidelity of detection is discussed. In Section El 
we briefly resume the results obtained and outline the prospects. 



2. Quadrature distributions of SPACS 



Generation of SPACS is due to injection of a coherent state \a) into the signal mode of an 
optical parametric amplifier. The stimulated emission of a single down-converted photon 
into the signal mode results in SPACS generation, which is trigged by the detection of 
a single photon in the idler mode of the amplifier. A time-domain balanced homodyne 
detector is then used to acquire quadrature data (see, e.g., the review [8]). 

The balanced homodyne detection is known to give access to quadratures Xg = 
Q cos# + Psin9, where [Q, P] = i and 9 £ [0, 2it] is a phase of an intense coherent light 
(the local oscillator). Once 6 is fixed, the distribution of quadratures is given by the 
optical tomogram w t h{X, 9) = (Xg\p\Xg), where p is the density operator of quantum 
state and X e \X e ) = X\X e ). 

Let p be a density operator of SPACS, then the tomogram w t h(X, 9) = (Xg\p\Xg) 
is easy to compute. However, it turns out that the experimentally measured quadrature 
distributions are smoother than the predicted ones and can differ significantly from 
them. This takes place due to losses and overall efficiency of detection rj < 1. One can 
make allowance for losses by introducing a fictitious beamsplitter with transmittivity rj 
in front of the ideal photodetectors (with sensitivity of 100%). Such an attenuation 
of the signal results in the following convolution relation between the quadrature 
distributions [30] : 

1 



w th (X, 9; rj) 



w th (Y,9)exp 



V 



V 



(y - ^xy 



dY. 



(1) 



vMl-r?) 

In the Schrodinger picture, the distribution w t h(X, 9; rf) is nothing else but 
the optical tomogram (XglS^pWXg) of the transformed state £ v [p], where E n [*] — 
YlT=o Ak(v) • A-liv) is a completely positive trace preserving map with the following 
operator-sum representation: Ak(i]) = Xlm=o ?7 m (l ~ v) k \m)(m + k\ (see, 

e.g., ISIED). 

Using formula ([TJ, one can calculate in explicit form the optical homodyne 
tomogram of a SPACS. Some algebra yields 

1 " r))[l + 4r)\a\ 2 sin 2 (9 -<p) 



w th (X, 9; a, rj) 



+ 2r] 



x exp 
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An analogue of tomogram (j2J) was first derived in the paper [15], where the authors 
used a slightly different commutation relation [Q,P] = |. The deduced tomogram ([2]) 
comprises two parameters: a = \a\e lip determines the coherent state |a) to which a single 
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photon is added, r\ is the overall efficiency of homodyne detection and characterizes the 
imperfection of measurement device. The overall efficiency includes transmission losses, 
mode matching and the intrinsic quantum efficiency of detectors. 



3. Estimation of parameters 



Our goal is to compare w t h{X, 9\ a, rj) with the experimentally measured distributions 
w ex (X, 9) and find parameters a = |a|e 1?i and rj resulting in the best fitting. In this sense, 
we perform a minimal distance estimation of the state parameter a and the detector 
parameter 77. In order to give this procedure more rigorous formulation with clearer 
physical meaning, we need to choose such a distance between distributions Wth(X, 9; a, rj) 
and w ex (X, 9) that is related with some fair distance between states p t h = £»?[p] an d p cx 
(satisfying metric requirements). Moreover, we are interested in such a distance between 
the states that could be operationally calculated via optical tomograms. Some aspects 
of appropriate distances were discussed in the paper [32]. The Hilbert-Schmidt distance 
D = A/Tr(p th — Pcx) 2 turns out to be suitable because it can be given by the following 
expression in terms of tomograms: 

1 r+oo r r+oo 

D 2 (a,r]) = - drrll dXdY cos[(X + Y)r\ 

7T J J J -00 



x 



d9 



w cx (X, 9) - w t h{X, 9; a, 77) 



w. 



-Y,9)-w th (-Y,9-a,r ] ) , (3) 



which can be readily deduced with the help of a formula for Tr pip 2 obtained in Ref . [33J . 
Similarly, the experimental error is evaluated by a slight modification of formulas in the 
paper [28], namely, 

1 r+00 p p+00 
A(D 2 ) = — drrll dXdY ^os[(X + Y)r] 

w ex {X, 9)w cx {-Y, 9) - w cx {X, 9 + tt)w cx {-Y, 9 + tt) 
+ 2w th (X,9)(w cx (Y,9 + ir) 



x 



d9 



-Y, 



(4) 



Formula (jlj) is based on the fact that the fair quadrature distributions satisfy the 
symmetry relation w(X, 9 + tt) = w(—X, 9). Experimentally measured distributions do 
not satisfy precisely this relation, and this gives rise to the error (jlj) which includes both 
systematic and statistical components (see details in the paper [28]). 



3.1. Results 

In this subsection, we estimate parameters |a|, <p, and rj for a particular set of 
experimental quadrature distributions. Phases of the local oscillator take discrete values 
{9j}^Li- For each fixed phase, the quadrature distribution is a histogram of 5321 values, 
with the bin width being chosen to guarantee the statistical confidence and prevent 
the data from undersampling |28j. Examples of experimental histograms are depicted 
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Figure 1. Experimental histograms w ex (X, 9) (blue discontinuous lines) and the 
closest theoretical quadrature distributions w t h{X,9) (red solid lines) of a SPACS for 
different phases: 9 = (a), 9 = 1.36 (b), 9 = 2.49 (c). 



in Fig. [TJ Thus, the data are presented in discrete form so the integrals in formulas 
(J3]) and (J4]) are calculated approximately by the trapezoid method [35J. The error of 
calculation is estimated in the paper [28] and is usually less than the experimental 
quantity (j4]). 

In our particular case, the minimization of the square of distance D 2 (\a\,ip,r]) 
results in D 2 = 0.0436 which is achieved at |a| opt = 0.81, ip opt = 3.14, 77 opt = 0.58. 
On substituting these parameters in formula ([2]), we can depict the closest theoretical 
quadrature distributions (see Fig. [1]). 

In order to evaluate the errors of estimated parameters |a| opt , y? op t, and i] opt we 
consider three cuts of the function D 2 (\a\,(p, rf) that cross at the point (|a| opt , <p pt, f] opt ). 
The values of function and their errors are shown in Fig. [2j Further, the error of an 
optimal parameter g op t can be evaluated as Ag/SNR, where Aq is the width of the 
corresponding function cut and SNR = (max/) 2 — minD 2 )/ max A(D 2 ) plays the role 
of signal to noise ratio. The errors evaluated in such a way give rise to the following 
results: |a| op t = 0.81 ± 0.03, ip opt = 3.14 ± 0.25, r] opt = 0.58 ± 0.02. The least precise 
parameter is the phase <p and this can be attributed to the relatively small mean number 
of photons (n) < 1 and imprecise fixing of the local oscillator phase 9. Improving control 
of this parameter would result in higher precision of parameters under estimation. 

4. Fidelity of detection 

Sometimes, the Hilbert-Schmidt distance between the states is not very representative 
because it can grow under the action of quantum operations (not monotone metric). In 
this case one exploits some other quantities, e.g., the Bures distance D B = a/2(1 — F), 
where F = Tr */ ^/pth"p cx y/fhh is Uhlmann's fidelity (see, e.g., the book [36])- The fidelity 
is difficult to express in operational way through quadrature distributions. Nevertheless, 
we can use recently found bounds for fidelity: the sub-fidelity E and the super-fidelity 
G satisfying E < F 2 < G and given by formulas [26l 1371 [38] 

E(pth, Pox) = Tr p th p cx + A/2[(Trp th p ex ) 2 - Tr p th Pe X PthPe X ], (5) 

G(p t h, Pcx) = Tr p th p cx + yj(l- Trp t 2 h ) (1 - Trp 2 x ). (6) 
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Figure 2. Square of the Hilbert-Schmidt distance vs. state and detector parameters 
in the vicinity of the global minimum: D 2 (\a\, ip opt , rj opt ) (a), D 2 (|a| op t, <p, rfopt) (b), 
D 2 (\a\opt,<Po P t,v) (c). 



The overlap Tr p t hp ex and purity Tr p 2 x are readily expressed through tomograms 
(see, e.g., [33]). It is worth noting that the purity can also be estimated by exploiting the 
covariant uncertainty relation |34J . As far as 4-product Tr PthPexPthPex is concerned, we 



can approximate it by Trp^ h . In fact, we have | Tr p th p cx p th p cx -Tr p^J < | Tr p cx p th p cx - 
Tr p^J < | Tr p 2 x — Tr p 2 J and can modify the sub-fidelity as follows: 



E'(p th , Pcx ) = Tr p th p ex + y / 2[(Tr p thPcx ) 2 - Trp* h - | Trp 2 x - Trp t 2 J]. (7) 

In order to be able to calculate the modified sub-fidelity ([7j) for SPACS, we find the 
following theoretical values: 



Tr Pt 2 h = !" 



2/7(1 - rj) 



Tr Pth 



4r;(l — 77) 2r] 2 (l-r]f 



(1 + M 2 ) 2 ' rtft (1 + |a| 2 ) 2 (1 + |a| 2 ) 4 ' 

Returning to the example considered earlier, we substitute the experimental data 
and the optimal theoretical values |a| op t — 0.81, <p opt = 3.14, 77 opt = 0.58 in formula 
OH]) and obtain the upper bound G(pth)Pex) = 0.98 ± 0.02. In our case, the direct 
calculation of sub-fidelity fl7]) turns out to be problematic because the confidence 
interval of the radicand is [-0.07; 0.05] (cf. 2[(Trp t 2 h ) 2 - Trp^] = 0.032). Thus, the 
calculation of square root is worthless. Therefore, the use of formula (J7|) is possible 
only with the data of very high precision (errors should be substantially less than 1%). 
Whenever this does not happen, one can use another lower bound E" = Tr pthPex < F 2 
(see, e.g., [26]). This lower bound is easy to calculate and in our case it equals 
E" = 0.81 ± 0.02. Consequently, the fidelity of our interest is bounded by the two- 
sided inequality 0.81 ± 0.02 < F 2 < 0.98 ± 0.02. 



5. Conclusions 



In order to estimate parameters of some prepared SPACS we developed the operational 
method whose essence was the comparison of experimental histograms with theoretically 
predicted quadrature distributions. The explicit form of theoretical distributions took 
into account the losses presented, which allowed us to infer not only the state parameter 
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a but also the parameter 77 describing the overall efficiency of homodyne detection. 
We discussed some practical issues concerning the easiest way to calculate the Hilbert- 
Schmidt distance and evaluate the errors of estimated parameters. The phase of the 
state turned out to be the least precise parameter, which could be ascribed to the small 
intensity of the signal mode and the errors in control of the local oscillator phase. Then 
we considered some operational techniques to determine the lower and upper bounds 
for fidelity of detection. We showed that, in practice, some of these bounds can be 
calculated only with highly precise data. 

The outlook for further research is to use high sensitivity of homodyne detection to 
trace all the stages of quantum state's life: its preparation, transformation via a quantum 
channel, and detection. Using appropriate theoretical models of these processes, one 
can determine the corresponding parameters. For instance, dark counts in the trigger 
detector result in mixing of the SPACS with a residual coherent state. In this case, 
the measured tomogram reads (1 — p)wspacs + V "^coherent , where p is a fraction of dark 
counts. The parameter p can be estimated by the same algorithm of comparing w ex and 

^th- 
in general, optical tomograms can be valuable information sources on equal footing 
with other state descriptions [39J. Improving the accuracy of homodyne detection, 
one can check the validity of more complicated quantum theories and observe new 
phenomena (see, e.g., [29]). The role of SPACS states for new experiments can be also 
dramatic because of their ability to exhibit properties ranging from classical to quantum 
ones for different intensities |14j . 
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